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Abstract 

By using the generalized Abel-Plana formula, we derive a summation formula for the 
series over the zeros of a combination of the associated Legendre functions with respect to 
the degree. The summation formula for the series over the zeros of the combination of the 
Bessel functions, previously discussed in the literature, is obtained as a limiting case. As 
an application we evaluate the Wightman function for a scalar field with general curva- 
ture coupling parameter in the region between concentric spherical shells on background of 
constant negative curvature space. For the Dirichlet boundary conditions the correspond- 
ing mode-sum contains series over the zeros of the combination of the associated Legendre 
functions. The application of the summation formula allows us to present the Wightman 
function in the form of the sum of two integrals. The first one corresponds to the Wightman 
function for the geometry of a single spherical shell and the second one is induced by the 
presence of the second shell. The boundary-induced part in the vacuum expectation value 
of the field squared is investigated. For points away from the boundaries the corresponding 
renormalization procedure is reduced to that for the boundary- free part. 

PACS numbers: 02.30.Gp, 03.70.-hk, 04.62.-hv 



1 Introduction 

The associated Legendre functions are an important class of special functions that appear in a 
wide range of problems of mathematical physics. The physical importance of these functions is 
related to the fact that they appear as solutions of the field theory equations in various situations. 
In particular, the radial parts of the solutions for the scalar, fermionic and electromagnetic 
wave equations on background of constant curvature spacetimes are expressed in terms of the 
associated Legendre functions (see, for instance, [B El [3]). The eigenfunctions in braneworld 
models with de Sitter and anti-de Sitter branes are also expressed in terms of these functions 
(see [1]). Motivated by this, in [5], by making use of the generalized Abel-Plana formula, we have 
derived a summation formula for the series over the zeros of the associated Legendre function 
of the first kind with respect to the degree (for the generalized Abel-Plana formula and its 
applications to physical problems see [H [71 [8] ) . This type of series is contained in the mode- 
sum for two-point functions of a quantum scalar field in background of a constant curvature 
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space with spherical boundary, on which the field obeys the Dirichlet boundary condition. The 
application of the summation formula allowed us to extract from the vacuum expectation values 
the part corresponding to the situation without boundary and to present the boundary-induced 
part in terms of rapidly convergent integral. 

In the corresponding problem with two concentric spherical boundaries, in the region between 
two spheres the eigenfunctions are the combination of the associated Legendre functions of the 
first and second kinds. The eigenfrequences are determined by the location of the zeros of this 
combination with respect to the degree. In the present paper, by specifying the functions in the 
generalized Abel-Plana formula, we obtain a summation formula for the series over these zeros. 
As in the case of the other Abel-Plana-type formulae, previously considered in the literature, 
this formula presents the sum of the series over the zeros of the combination of the associated 
Legendre function in the form of the sum of two integrals. In boundary-value problems with 
two boundaries the first integral corresponds to the situation when one of the boundaries is 
absent and the second one presents the part induced by the second boundary. For a large 
class of functions the latter is rapidly convergent and, in particular, is useful for the numerical 
evaluations of the corresponding physical characteristics. 

The paper is organized as follows. In section [21 by specifying the functions in the generalized 
Abel-Plana formula we derive a formula for the summation of series over zeros of the combination 
of the associated Legendre functions with respect to the degree. In section [3l special cases of 
this summation formula are considered. First, as a partial check we show that as a special case 
the standard Abel-Plana formula is obtained. Then we show that from the summation formula 
discussed in section [21 as a limiting case the formula is obtained for the summation of the series 
over the zeros of the combinations of the Bessel functions, previously derived in [6] . A physical 
application is given in section [4l where the positive frequency Wightman function for a scalar 
field is evaluated in the region between two spherical boundaries on background of a negative 
constant curvature space. It is assumed that the field obeys Dirichlet boundary condition on 
the spherical shells. The use of the summation formula from section [21 allows us to extract from 
the vacuum expectation value the part corresponding to the geometry where the outer sphere 
is absent. The part induced by the latter is presented in terms of an integral, which is rapidly 
convergent in the coincidence limit for points away from the sphere. The main results of the 
paper are summarized in section \5\ In appendix [A] the formula for the normalization integral 
is derived and we show that the zeros of the combination of the associated Legendre functions 
with respect to the degree are simple. 



2 Summation formula 

Let z = Zk, k = 1,2, . . ., he zeros of the function 
in the right-half plane of the complex variable z: 

K,M = ^- (2) 

In ([T|), associated Legendre function of the first kind (in this paper the definition 

of the associated Legendre functions follows that given in [9]). In the discussion below we will 
assume that u,v > 1. The expression in the numerator of ([T]) has simple zeros for integer values 
of fj, and the function X^^{u, v) is regular at these points. Since one has the property Xi7^(u, v) = 
Xy{u, v), without loss of generality, we consider the parameter ^ being non-negative, /i ^ 0. For 
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given values u, v, and /i the function X^^{u, v) has an infinity of real zeros. From the asymptotic 
formula for the associated Legendre functions we can see that for z +00 one has 



2sm[{r]y - r]u) z] 



7rz\/sinh r]u sinh rjy ' 
where r]u and T]y are defined as 

U = cosh 7]u, V = cosh 7]y . 

From here we obtain the asymptotic expression for large zeros: 

Zk ^ irk/ {r]y - r]u) . 



(3) 



(4) 



(5) 



In general, the zeros are functions of the parameters n, v, and fi: = Zk{u,v, fi). By 
taking into account that for the associated Legendre function one has P^^_ij2{'^) = ^^-1/2^''^)^ 
we see that X^^{u, v) = Xi}{u, v). Hence, the points z = —z^ are zeros of the function X^^{u^ v) 
as well. In Appendix [X] we show that the zeros z = z^ are simple and under the conditions 
specified above the function X|i,(n, has no zeros which are not real. We will assume that 
Zk are arranged in ascending order of magnitude. Note that the function Xy{u,v) can also be 
expressed in terms of the combination 



as 



Xl^{u,v) 



2 r(iy-/u + l/2) 

Tj-(,ifi7T r(i/ + /u + 1/2) 



(6) 



(7) 



where Q'^_i^2(^) associated Legendre function of the second kind and r(x) is the gamma 

function. 

A summation formula for the series over Zk can be derived by using the generalized Abel- 
Plana formula [6j (see also [ZllS])- For functions f{z) and g{z) meromorphic in the strip a ^ x 
of the complex plane z = x + iy this formula has the form 



lim 

6— >oo 



p.v. 



dx f{x) 



^ Res g{z) 

^ 2— Zg J. 



TTl 



Ea{zf^k) Res f{z) 



a+ioo 



dz [giz)+a{z)fiz)] 



(8) 



where cr{z) = sgn(Imz) and p.v. means the principal value of the integral. In this formula, Zf^^ 
and Zg^k are the positions of the poles of the functions f{z) and g{z) in the strip a < x < b. As 
functions f{z) and 17(2;) in formula ([8]) we choose 



fiz) 
9{z) 



h{z) 



r{iz + fi+ 1/2) vr^e^^^'^i sinh( 



ZTT 



^Qiz-i/2(^)Q-iz-i/2^^) r (iz - ^ + 1/2) cos[(iz - /u)7r] 



Qli,^l/2i^) 



+ 



Qiz~l/2 



(n) 



hiz) 



2Xf,(n,. 



(9) 



where h{z) is a meromorphic function for a ^ Re z ^ b. The combinations appearing on the 
left-hand side of formula ([8]) are presented in the form 



9{z)±f{z) 



Qt^^-i/2(^) h{z) 



(10) 
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Note that the function g{z) has simple poles at the zeros of the function ([T]). 

The conditions for the generalized Abel-Plana formula ([8]), formulated in terms of the func- 
tion h{z), take the form 



lim 

ui— »oo 



lim 



, QTiz-iM"") h{z) 

a±iw 



dz 



{u)X>^^{u,v) 
(v) h{z) 



0, 



0. 



(11) 



With the help of the asymptotic formulae for the associated Legendre functions, we can see that 
these conditions are satisfied if the function h{z) is restricted by the constraint 



\h{z)\ < x-2Me(a;)e^('?--'?")y, z = x + iy, \z\ 



oo, 



(12) 



uniformly in any finite interval of x, where c < 2, e(x) — > for x — > +oo. 

Now, after the substitution of the functions ([9|) into formula ([8]), we see that for a function 
h{z) meromorphic in the half-plane Rez ^ a and satisfying condition ()12p . the following formula 
takes place 



lim 



E 



h(z) e£-i/2(^') 



Z=Zk 



■"+*~ Q-a(z)^z-l/2i^) h{z) 

Q^. ...... .MXf 



H — p.v. 

TT 



dx f{x) + r[h{z)] 



(13) 

-(T{z)iz-l/2^ 

where the function f{z) is defined by the relation ([9]). In this formula we have introduced the 
notation 

r[h{z)] = ^ Res 



Z = Zh.k 



fc,Im2h t=0 



Res 



Q-a{zk)iz-l/2^'"^ h{z) 
Q-aiz,)^z~l/2(^) ^izi^^^) 



^Hz-1 



h{z) 



12^ 



^r.(^,-),t^QL-i/2(") 



(14) 



with z/i^fc being the positions of the poles for the function h{z\ On the left-hand side of ()13p . 
one has z^-i < a < Zm, Zn < b < Zn+i and in (fH|) the summation goes over the poles Zh^k in 
the strip a < Re z < b. Note that one has the relations 



^7, 



iz-l/2y 



-^12-1/2 



^1.2 



iz-l/2^ 



jz-1/2 



(^) 



^:^/2(-) 



^;^/2(-) 



(15) 



and in the summation of the first term in figure braces of (jl3p we can replace the ratio of the 
associated Legendre functions of the second kind by the ratio of the functions of the first kind. 

A useful form of the summation formula (|13p is obtained in the limit a — > 0. In this limit, we 
see that for a function h[z) meromorphic in the half-plane Rez ^ and satisfying the condition 
(|12p the following formula holds 



h{z) 



42-1/2 



(v) 



h(x) 



u 



vre 



2i/i7r 



-p.v. 



dx 



Z = Zk 



^ix-1/2 



(")Q-ix-l/2(^) 



r {ix + ^ + 1/2) sinh(x7r) 
r {ix — ^ + 1/2) cos[(ix — /i)7r] 

°° , Qx-i/2^'") h{xe^'/^) + /i(xe-^*/2) 



dx 



(16) 
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For large values x ^ 1, for the associated Legendre functions we have 



By using these formulae and the relation d?]), for the corresponding asymptotic behavior of the 
function Xx{u,v) one finds 

JVv-Vn)x 

Xliiu,v)f^ . . ■ . , . (18) 

TTX V smh r]u smh rj^ 

Prom these asymptotic formulae it follows that under the condition ()12p for the function h{z), 
the second integral on the right-hand side of formula (|16p exponentially converges in the upper 
limit. 

If the function h{z) has poles on the positive real axis, it is assumed that the first integral 
on the right-hand side converges in the sense of the principal value. From the derivation of 
(|16|) it follows that this formula may be extended to the case of some functions h{z) having 
branch-points on the imaginary axis, for example, having the form h(z) = hi{z)/{z^ + c^)^/^, 
where hi{z) is a meromorphic function. This type of function appears in the physical example 
discussed in sectional Special cases of formula (fT6|) are considered in the next section. 

Another generalization of formula (jl6p can be given for a class of functions h{z) having purely 
imaginary poles at the points z = ±iyk, Uk > 0, k = 1,2, .. ., and at the origin z = yo = 0. We 
assume that the function h{z) satisfies the condition 

h{z) = -h{ze^'^') + o{{z - cjfc)"^), z ak, Ok = 0, iyk- (19) 

Let us denote by Cp{ak) the right half of the circle with radius p and with the center at the 
point cTfc, described in the positive direction. Similarly, we denote by 7^ and 7" the upper and 
lower halves of the semicircle in the right half-plane with radius p and with the center at the 
point 2: = 0, described in the positive direction with respect to this point. Now, in the limit 
a ^ the right-hand side of (jlSh can be presented in the form 



plus the sum of the integrals along the straight segments (±i(yfc_i -|- p),±i{yk — p)) of the 
imaginary axis between the poles. In the limit p ^ the sum of the integrals along the straight 
segments of the imaginary axis gives the principal value of the last integral on the right-hand 
side of (jl6[) . In the terms of ()20p with a = — we introduce a new integration variable z' = ze" 
By using the relation (jl9p . the expression (j20p is presented in the form 



(l-5oa,/2)Res 



Q-i,-l/2i^) hjz) 



(21) 



plus the part which vanishes in the limit /? — > 0. As a result, formula (|16p is extended for 
functions having purely imaginary poles and satisfying condition ()19p . For this, on the right- 
hand side of ([16]) we have to add the sum of residues ([2T]) at these poles and take the principal 
value of the second integral on the right-hand side. The latter exists due to condition (|19p . 
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3 Special cases 



First we consider the case fi = 1/2. For the corresponding associated Legendre functions one 
has 

>-i/2,_^^, ../l^inM^^ pi/2 
vr z\ 



P- ^/.(coshr,) = V-;^^, P^.y^i^oshrj) 
By making use of these formulae we find 



2 cosh(zr7) 
TT y^smh?] 



(22) 



2 sm[z{r]^ - r/„)] 
vr z-^ sinh sinh rj^ 



(23) 



Hence, in this case for the zeros Zj^ one has zj^ — irkfi^Tjy — 7]u). Introducing a new function F{z) 
in accordance with the relation zh{z) = F{z{r]y — r/„)/7r), from the formula (|16p we obtain the 
Abel-Plana summation formula in its standard form: 



k=l 



dx F{x) + i 



°° , F(ix 
ax 



F{-ix) 



2nx 



1 



(24) 



where the first term on the right-hand side comes from the residue term at fifc = in (|2ip. 

Now let us show that from formula (jl6p . as a special case, a summation formula is obtained 
for the series over zeros of the combination of cylinder functions. First of all, by making use of 
formulae 



hmJszf''P^/;_^^^icosh{X/s)) = J±^{\2 



(25) 



with J^{r}) being the Bessel function of the first kind, we can see that the following relation 
holds: 

lim Xj'^^(cosh(Au/s),cosh(A^/s)) = C^(AuZ, A^z), (26) 



where 



C^{\uZ^\vz) — J^{Xuz)Yfj_(Xvz) — Yfj_[Xuz)J^[Xyz). 



(27) 



Note that, instead of the function J-^[z) we have introduced the Neumann function Y^{z). 
Hence, in the limit s ^ oo from (jl6p we obtain the summation formula for the series over zeros 
z = Xf^^k, k = 1,2,..., of the function C^{XuZ, X^z). For this, first we rewrite formula (fT6]) 
making the replacements z — > sz, x — > sx, in both sides of this formula including the terms in 
r[/i(z)], and we take u = cosh(Au/s), v = cosh(Ai,/s). Introducing a new function F{z) = h{sz), 
in the limit s — > +oo we find the formula 



E 



F{z) MX^z) 



^ dzC^{XuZ, Xyz) Jfj,{Xq^ 



1 



fj.,k 



vr 



-p.v 



dx ■ 



F{x) 



-Tc[F{z)] - - 



Jl{Xux) +Y^{Xux) 
F{xe^'/'^)+F{xe-"'/'^) 



K^{Xux) K^{Xux)Ifj,{Xyx) - I^{Xux)K^{Xvx) 



where I^{x) and K^{x) are the modified Bessel functions and 



rc[F{z) 



vr > Res 

^-^ Tm y „ , - 



J^,{Xuz)J^{X^z) + Y^{Xuz)Y^{X^z) F{z) 



mXux) +Y^{Xux) 



C^{XuZ, Xyz) 



+vr "S^ Res 



Hi-\x,z) 



F{z) 



hP{Xuz)C^{XuZ,Xvz) 



(28) 



(29) 
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In deriving (j28|) we have also used the formulae 



and the relation 



lim i^^^P;^^ {cosh{x/u)) = I±Jx), (30) 

I/— >oo 

lim u~^'Q'^{cosh{x/u)) = e*'^^K„(x), 



Note that from (1301) it follows that 



, z = X^^k- (31) 



2 

lim X^3,(cosh(A„/s),cosh(At,/s)) = - [Ka{Xux)Ia{Xvx) - iJXux)Ka(Xvx)] . (32) 

s—>-+oo TT 

Formula (p8]) is a special case of the result derived in ^ (see also, [8]). Physical applications of 
this formula are given in \10\ [TT] . 

4 Vacuum polarization by concentric spherical boundaries in a 
constant curvature space 



In this section we give a physical application of the summation formula (jl6p . Consider a scalar 
field '^{x) on background of the space with constant negative curvature described by the line 
element 

ds'^ = dt^ - [dr"^ + sinh^ ridd"^ + sin^ Odcf)^)] , (33) 
where a is a constant. The field equation has the form 

ViV' + m2 + ^r) ^{x) = 0, (34) 



where M is the mass of the field quanta, is the curvature coupling parameter and for the Ricci 
scalar one has R = -6a-2. We wih assume that the field operator satisfies Dirichlet boundary 
conditions on two concentric spherical shells with radii r = ri and r = r2, ri < r2, 

97(x)|r=ri,2 = 0. (35) 

The boundary conditions modify the spectrum of the zero-point fiuctuations and, as a result 
of this modification, the physical properties of the vacuum are changed. Among the most 
important characteristics of these properties are the expectation values of quantities bilinear in 
the field operator such as the field squared and the energy-momentum tensor. These expectation 
values are obtained from two-point functions in the coincidence limit of the arguments. As a 
two-point function here we will consider the positive frequency Wightman function. Other two- 
point functions are evaluated in a similar way. Expanding the field operator over the complete 
set {ipa{x), 93* (j;)} of classical solutions to the field equation satisfying the boundary conditions 
(I35p . the Wightman function is presented in the form of the following mode-sum 

Wix, x') = {0\ip{xMx')\0) = ^a{x)^l{x'), (36) 

a 

where |0) is the amplitude of the vacuum state and a is a set of quantum numbers specifying 
the solution. 
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In accordance with the spherical symmetry of the problem under consideration, the eigen- 
functions for the scalar field can be presented in the factorized form 

^o.{x) = Z{r)Yi^{9,<t))e-'^\ (37) 

where Yiy^iO, (p) are the spherical harmonics with / = 0, 1, 2, . . ., —l^m^l. The equation for 
the radial function is obtained from the field equation (|34|) and has the form 



1 d f 2 dZ \ 

smh r— — + 



sinh^ rdr \ dr 



sinh^ r 



Z = 0. (38) 



In the region between the spherical shells the solution of equation (p8]) is expressed in terms of 
the associated Legendre function as 

„-l-l/2, ^ d'+1/2 / X 

^^^^ ^ Cl-P».-l/2 (U) + C2P,,_i/2{u) 

\J sinh r 

with integration constants c\ and C2 and the notations 

= {u? - M'^)a^ u = cosh r. (39) 

From the boundary condition on the inner sphere we find 



and, hence, 



, Ui = coshrj, i = 1,2, (40) 



Zir) = C ^-'P^"^ . (41) 
Vsmhr 

where Ca is the normalization constant and the function X-^^''^(iii, u) is defined by ([T]). From 
the boundary condition on the outer sphere we see that the eigenvalues for z are solutions of 
the equation 

Xlt'^\u,,u2) = 0. (42) 
As a result, the eigenfunctions have the form 

^Ux) = -^xl^'/\uuu)YUe,<t>)e-'^', (43) 
V smh r 

and, hence, z = z^, A; = 1, 2, . . ., in the notations of section[2j The corresponding eigenfrequencies 
are related to these zeros by the formula 

u;l = ojHzk) = {zl + l-60/a^ + M^. (44) 

Hence, the set a of the quantum numbers is specified to a = (/, m, k). 

The coefficient Ca in (|43|) is determined from the orthonormalization condition for the eigen- 
functions: ^ 

d^x y^\ipa{x)ip*^,{x) = (45) 

where the integration goes over the region between the spherical shells. Making use of the 
integration formula given in Appendix |A] and the boundary conditions, for this coefficient we 
find 

Ca' = a'^^iul - l)[d.x':''\u,,u,)]duX'^''\u,,u), (46) 
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with z = Zk, u = U2- By using the Wronskian relation for the associated Legendre functions, 
it can be seen that 

p/ + l/2 / N 

[^«^-. (ni,n)]„ = - 2_. m/2 , (48) 

Upon the substitution this into (146 p . the normahzation coefficient is written in the equivalent 
form 

p/ + l/2 / N 

1/ 



^-(2(^2) 



(49) 

Note that the ratio of the gamma functions in this formula can also be presented in the form 

r(z^, + / + i) ^ 1 _ ^ _ i/2)]|r(izfe + 1 + (50) 

Substituting the eigenfunctions into the mode-sum formula (|36p and using the addition 
theorem for the spherical harmonics, for the Wightman function one finds 



^ 00 



(2; + l)Pf(cos7) 
fl'^ V sinh r sinh r' 

00 



(51) 



where At = t — t' and n' = coshr'. In (j5ip . P/(cos7) is the Legendre polynomial and 

cos 7 = cos cos 9' + sin sin 6' cos((/) — (/)')• (52) 

As the expressions for the zeros z^ are not explicitly known, formula (j5ip for the Wightman 
function is not convenient. In addition, the terms in the sum are highly oscillatory for large 
values of quantum numbers. 

For the further evaluation of the Wightman function we apply to the series over k the 
summation formula (jl6p with u = ui and v = U2, taking in this formula 

h{z) = zXl'^^''^{ui,u)xl'^^''^{ui,u) — — , (53) 

Ll![Z) 

where the function uj{z) is defined by (j44p . The function (jSSp has no poles in the right-half plane 
and, hence, r[/i(z)] = 0. The corresponding conditions are satisfied if r + r' + At/a < 2r2. In 
particular, this is the case in the coincidence limit t = t' for the region under consideration. For 
the function (I53p the part of the integral on the right-hand side of formula (I16p over the region 
(0, xm) vanishes, and for the Wightman function one finds 



W{x,x') = Wi{x,x') -\ ^ ' ) " 

Jl^ V sinh r sinh r' Jx 



00 

dx X 



CV2(^2) xi+^/^(ni,n)xi+^/^(ni,nO ^oMyJ - ^lA^/a) 
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where we have defined 

XM = ^jM'^a^ + 1 - (55) 
In formula (|54p . the first term on the right-hand side is given by 

Wxix.x) = ---"-^ V^^^^i4==T=^ / c^x X sinh(x7r) 
o^o vsmhrsmhr' Jq 



This function does not depend on the outer sphere radius whereas the second term in ([5i 
vanishes in the limit r2 — > oo. Hence, the two-point function given by ()56p is the Wightman 
function for a scalar field in background spacetime described by the line element (j33p outside a 
single sphere with radius ri on which the field obeys Dirichlet boundary condition. This can also 
be seen by the direct evaluation using the corresponding eigenfunctions. Thus, we can interpret 
the second term on the right-hand side of (|54p as the part in the Wightman function induced 
by the presence of the outer sphere. 

An alternative form for the function (j56p is obtained by making use of the identity 

X .^ (m,M)X.^ (m,n) _ _^p'i^\i2, . -1-XI2, 



'3ix-l/2(^l)Q-*x-l/2(^l) 



TT' 



4« p-«-l/2. ^ Qaix-l/2NQa.x-l/2(^) , . 

Z^ix-l/2 ^-i-1/2 , . • 

±1 Q<x*x- 72(^1) 



Substituting (j57p into (j56p . we can see that the part with the first term on the right of for- 
mula (15] 



..r r 1 ^ (2Z + l)P,(cos7) /"^ , , 

Wo[x,x) = — / — 1^^^=^^^^ / dx X smh[7rx) 

Vsinhrsinhr' Jo 

x\T{ix + / + 1)|2f;^Y/2 (^oshr)FiY/2 (^oshr')^^, (58) 

is the Wightman function for a scalar field on background of the constant curvature space 
without boundaries (see [5])- In the part with the second term on the right-hand side of formula 
(j57p we rotate the contour of integration over x by the angle tt/2 for the term with a = —1 and 
by the angle — 7r/2 for the term with a = 1. As a result, the exterior Wightman function for a 
single spherical boundary is presented in the decomposed form 

/X /N i V^/ ,,(2/ + l)P/(cos7) /""^ , r(3; + / + l) 

W,{x,x') = Wo{x,x') - -^Y.^-iy-T^ u u / dxx^ 

47r^a^ ^ vsmhrsmhr' Jxm T{x-l) 

P^-ij'i'^ i'^i) -1-1/2 -1-1/2 cosh{J x^ — x\jAt/a) 
X Qx-1/2 (cosh r)Q^_^// (cosh/) \ , (59) 

where the second term on the right-hand side is induced by the spherical boundary. The Wight- 
man function for the region inside a single spherical shell is investigated in [53- The corresponding 
expression is obtained from ([59|) by the replacements ^x-\/2 ™ second term on 

the right of this formula. 
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Taking the limit o ^ oo with fixed ar = R, from the formulae given above we obtain the 
corresponding results for spherical boundaries in the Minkowski spacetime with radii Ri = ari 
and i?2 = (^1^2- Note that in this limit one has xm = ciM and the result does not depend on the 
curvature coupling parameter. Introducing a new integration variable y = x/a and using the 
asymptotic formula for the gamma function for large values of the argument, from (j54p we find 



(2/ + l)Pi(cos7) cosh(^y2 _ M^/\t) 



Ki+i/2{R2y) Gi+y^jRiy. Ry)Gi+y2{Riy, R'y) , . 

Ki+,/2{Riy) Gi+y2{Riy,R2y) ' ^ ' 

where we have introduced the notation Gu{x,y) = Ky{x)Iy{y) — K^{y)I^{x). The first term 
on the right-hand side of formula (j60p is the Wightman function in the region outside a sin- 
gle spherical boundary with radius Ri in the Minkowski bulk. This function is given by the 
expression 



w(M)/ /N TI7{M)/ /N ^ {21 + 1) Pi {cos j) f 

W} '(x,x) = VFn '(x,x)-y , / 



oo 

dyy 

M 



i+i/2> yt I+1/2V i"k,+,/2(Ki!/) Vs" - 

Expressions (j60p and (j6ip are special cases of the general formulae given in [10] for a scalar field 
with Robin boundary conditions in arbitrary number of spatial dimensions. 

The vacuum expectation value of the field squared is obtained from the Wightman function 
taking the coincidence limit of the arguments. This limit is divergent and some renormalization 
procedure is necessary. Here the important point is that for points outside the spherical shells 
the local geometry is the same as for the case of without boundaries and, hence, the structure 
of the divergences is the same as well. This is also directly seen from formulae ()54p and (j59p . 
where the second terms on the right-hand sides are finite in the coincidence limit. Since in these 
formulae we have already explicitly subtracted the boundary- free part, the renormalization is 
reduced to that for the geometry without boundaries. In this way for the renormalized vacuum 
expectation value of the field squared one has 



ioa smnr j^^^ ^3.2 _ 3.2^ 

r(x + / + 1) QL-Y/'2(^2) u)]2 



(62) 



where the first term on the right-hand side is the corresponding quantity outside a spherical 
boundary with radius ri in the constant negative curvature space without boundaries and the 
second one is induced by the presence of the second spherical shell with the radius r2. For the 
first term one has 

{^)l,ren = (V' )o,ren " >^ ^353 -^-7— dx X 



^ Air'^a'^ sinhr 



-1-1/2, 12 



(63) 



11 



where (¥'^)o,ren is the vacuum expectation value for the field squared in the constant negative 
curvature space without boundaries and the second one is induced by the presence of a single 
spherical shell with radius ri. Note that the corresponding formula for the vacuum expectation 

^ 1/2 

value inside a spherical shell (see [5]) is obtained from (f63ll by the replacements -Pj._;^/2 ^ 

^x-\l2 second term on the right-hand side. 

The physical example discussed in this section demonstrates the advantages for the applica- 
tion of the Abel-Plana-type formulae in the evaluation of the expectation values of local physical 
observables in the presence of boundaries. For the summation of the corresponding mode-sums 
the explicit form of the eigenfrequencies is not necessary and the part corresponding to the 
boundary-free space is explicitly extracted. Further, the boundary induced parts are presented 
in the form of integrals which rapidly converge and are finite in the coincidence limit for points 
away from the boundaries. In this way the renormalization procedure for local physical observ- 
ables is reduced to that in quantum field theory without boundaries. Methods for the evaluation 
of global characteristics of the vacuum, such as the total Casimir energy, in problems where the 
eigenmodes are given implicitly as zeros of a given function, are described in references |12j . 



5 Conclusion 

The associated Legendre functions arise in many problems of mathematical physics. By making 
use of the generalized Abel-Plana formula, we have derived summation formula (jl6p for the series 
over the zeros of the combination ([T]) of the associated Legendre functions with respect to the 
degree. This formula is valid for functions meromorphic in the right half-plane and obeying 
condition (I12p . The summation formula may be extended to a class of functions having purely 
imaginary poles and satisfying the condition (|19p. For this, on the right-hand side of (|16p we have 
to add the sum of residues (j2ip and take the principal value of the second integral on the right- 
hand side. Using formula ()16p . the difference between the sum over the zeros of the combination 
of the associated Legendre functions and the corresponding integral is presented in terms of an 
integral involving the Legendre associated functions with real values of the degree plus residue 
terms. For a large class of functions h{z) this integral converges exponentially fast and, in 
particular, is useful for numerical calculations. The Abel-Plana summation formula is obtained 
as a special case of formula (|16p with = 1/2 and for an analytic function h{z). Applying the 
summation formula for the series over the zeros of the function Xj^(cosh(Au/s), cosh(A^/s)) and 
taking the limit s ^ oo, we have obtained formula (j28p for the summation of the series over 
zeros of the combination of the Bessel functions. The latter is a special case of the formula, 
previously derived in [6]. 

A physical application of the summation formula is given in section [H For a quantum scalar 
field with the general curvature parameter we have evaluated the positive frequency Wightman 
function and the vacuum expectation value of the field squared for the geometry of concentric 
spherical shells in a constant negative curvature space. The Dirichlet boundary conditions on 
both shells are assumed. In the region between the shells the eigenfunctions have the form ()43p 
and the corresponding eigenfrequencies are related to the zeros of the function x'^^^^(ui, 
by the formula (j44p . For the evaluation of the corresponding series in the mode-sum (jSip for the 
Wightman function we apply summation formula ()16p with the function given by ()53p . As 
a result this function is presented in the decomposed form ()54p , where the first term on the right 
is the Wightman function for the region outside a single spherical boundary and the second one 
is induced by the presence of the outer sphere. By making use of the identity (j57p . we have 
presented the single shell Wightman function as a sum of two terms, formula (j59p . The first one 
is the corresponding function in the constant curvature space without boundaries and the second 
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one is induced by the shell. For points away from the shell the latter is finite in the coincidence 
limit and can be directly used for the evaluation of the boundary-induced part in the vacuum 
expectation value of the field squared. The renormalization is necessary for the boundary-free 
part only and this procedure is the same as that in quantum field theory without boundaries. 
In the region between the spherical shells the vacuum expectation value of the field squared 
is presented in the form (|62p . where the first term on the right-hand side is the corresponding 
quantity outside a spherical boundary and is given by the expression (fMI) . 
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A On the zeros of the function f) 

In this appendix we show that the zeros z = are simple and real. First, we note that the 
functions -Pj^^^i/2(^) satisfy the Legendre equation with v = iz — 1/2 and, hence, the function 
Xf^{u,v) is a solution of the Legendre equation for the same value v with respect to both 
arguments. As a result, by making use of the differential equation for the associated Legendre 
functions it can be seen that the following integration formula takes place 

/ duX'^,{u,v)Xii{u,v) = [l-u )^ 2 • (64) 

Taking the limit v' ^ v and applying the I'Hopital's rule for the right-hand side, from this 
formula we find 

^yu\Xl{u,v)f = -^^{[5,Xf,(n,r;)]9„Xf,(n,^) -Xf,(^,T;)9.9.Xf,(n,^)}. (65) 

By taking into account the relation X'^-^{u,v) = X-^^{u,v), we see that for real z one has 
[Xj^('u)]^ = |Xj^(m)P and the integral on the left-hand side of ([65|) is positive. Now from (j65|) it 
follows that [dzX^^{u,v)]z=Zk / 0) and, hence, the zeros Zk are simple. 

Now let us show that all zeros of the function X^^{u,v) are real. Suppose that z = A is a 
zero of X-t,{u,v) which is not real. As the function X-^^{u,v) has no pure imaginary zeros, A is 
not a pure imaginary. If A* is the complex conjugate to A, then it is also a zero of X^^{u,v), 
because X^^,{u,v) = [X^^{u,v)]* . As a result, from formula (p4l) we find 

\vXl'^,{u,v)Xl'^{u,v)=0. (66) 

We have obtained a contradiction, since the integrand on the left hand-side is positive. Hence 
the number A cannot exist and the function X^^{u) has no zeros which are not real. 
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